Absence of anisotropic universal transport in YBCO 
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There exists significant in- plane anisotropy between a and b axis for various properties in YBCO. 
However recent thermal conductivity measurement by Chiao et al. which confirms previous mi- 
crowave conductivity measurement by Zhang et al, shows no obvious anisotropy in the context of 
universal transport. We give a possible explanation of why the anisotropy is seen in most properties 
but not seen in the universal transport. 
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PACS numbers: 78.30 .-j, 74.62.Dh, 74.25.Gz 

The topic of anisotropic superconductors is currently 
of great interest in relation to high-temperature super- 
conductivity. For example, in the most studied material 
YBCO, due to the orthorhombic nature of its band struc- 
ture, it exhibits strong in-plane anisotropy between some 
measurable quantities along the a and b axis. Anisotropic 
behavior is confirmed in the magnetic penetration depth 
jiormal-state resistivity p || , optical conductivity 
and thermal conductivity k || . Consistently a 
value of roughly two is found for the ratio of the carrier 
density to the effective electron mass (or the square of 
plasmon frequency) between the 6-axis and the a-axis. 
This implies, in regard to YBCO, that the contribution 
from the CuO chain is crucial and of the same order as 
the Cu02 plane contribution. 

It is well known, following a Drude-like model, that 
the normal-state DC electrical conductivity is given by 
<7„ = ne 2 r Jm with n the carrier density, m the effective 
mass of the carrier, and r the average scattering time. 
When the system undergoes a phase transition into a su- 
perconducting d x 2_ y 2 state, a universal microwave con- 
ductivity cr" = ne 2 /irmAg (Ao is the maximumgap on 
Fermi surface) is predicted at low temperatures |6| . This 
universality arises due to the cancellation between the 
finite impurity-induced quasiparticle density of states at 
zero energy and the scattering time. Universal features 
have been confirmed in YBCO by Zhang et al. [fjj for the 
microwave conductivity and by Taillefer et al. H] for the 
thermal conductivity. Clearly a n and a® share a com- 
mon dependence on n/m so that one expects a similar 
in-plane a-b axis anisotropy for them. Apart from the de- 
pendence on n/m, an anisotropic a n can also arise from 
an anisotropic scattering time r, while the universal a®, 
by the meaning of "universal" , is independent of the scat- 
tering rate. 

More recent thermal conductivity measurements of 
Chiao et al. Q which confirm previous microwave con- 
ductivity measurements of Zhang et al. B on YBCO have 
observed the predicted low-temperature universality || 
in both a and b directions, but observed no significant 
anisotropy between the two. This is in contradiction with 



the large anisotropy observed in many other quantities. 

In this paper, we propose a scenario that resolves the 
discrepancy mentioned above. In reality, YBCO has 
three conducting layers (two Cu02 plane and one CuO 
chain layer) within a unit cell. There is strong evidence 
|^| |l2| that the gap in all three layers is mainly d x 2_ y 2- 
wave like with some small admixture of an s-wave com- 
ponent |Il|,[l4]]. The sub-dominant component arises be- 
cause of the orthorhombic nature of YBCO. The Fermi 
surfaces are quite different in the different layers. For 
isolated layers, it is tetragonally two-dimensional for the 
two Cu02 planes, while it is highly one-dimensional for 
the CuO chain layer. 
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FIG. 1. Fermi surface for Cu02 plane (a) and CuO chain 
(b) layers in YBCO. Frames (c) and (d) are the d x 2 2-wave 
gaps seen on the Fermi surfaces for the plane and chain layers. 

A model that captures the physics in a simple way is 
the following. We assume all three layers in YBCO to 
have a circular Fermi surface as is usually done and as- 
sume a d x 2_ y 2-w&ve gap, described by Ak = Ao cos(20) 
as usual with <f> a polar angle in the 2D Brillouin zone. 
But to model for the single chain layer, we introduce 
a cutoff angle <fi c (see Fig. R1) which corresponds to an 
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open circular Fermi surface and simulates the highly one- 
dimensional nature of the chain band. The Fermi surface 
for chain band is shown as in Fig. |l|(b) with the CuO 
chain taken to be along the k v axis. A cutoff angle, 
4> c > t/4, corresponds to the angular position where a 
minimum gap starts to appear on the Fermi surface of the 
chain layer. The measured properties will correspond to 
the sum of the contributions from the two plane layers 
(associated with a complete Fermi-surface integration) 
and the one chain layer (associated with an incomplete 
Fermi-surface integration) . 

Measurements on YBCO can be divided into two cate- 
gories. The ones which are strongly dependent on seeing 
the gap nodes on Fermi surface (e.g., <r°) and another's 
which arc not too dependent such as the value of the 
zero-temperature penetration depth A(0). In the former 
case, due to the cutoff angle, the gap nodes may not cross 
anywhere on the entire Fermi surface for the chain band, 
and anisotropy is consequently suppressed because of the 
reduction of universal conductivity in the chain layer. 

Normal-State Transport. — Assuming that the scatter- 
ing rate is not dependent on the direction of the momen- 
tum k, the ratio of normal chain conductivity along the 
a and b axis is then given by (v is the Fermi velocity) 
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A, denotes an incomplete 



where, (A) c = (2/tt) / ( 
Fermi-surface average for the chain layer. In Fig. |2|, we 
plot the ratio, as a function of the cutoff angle <p c , given 
in Eq. (|l|). It is seen, as expected, that a h n > a® when an 
incomplete Fermi-surface integration is carried out be- 
cause (Vy) c has larger components than (vp c . When <fi c 
is large corresponding to a flatter band, u h n 3> a". For 
example in the case 4> c — 45°, cr^/cr^ ~ 4. This, in com- 
bination with the plane contribution, conforms with the 
large observed anisotropy between the the a and 6-axis 
resistivities ||. In a strictly ID model for the chains, a a 
would be zero. Here it is finite but becomes small when 
(j) c is large which approaches the chain case more closely 
and the model is good enough for qualitative arguments. 

Similar arguments can be applied to the anisotropy 
exhibited in the zero-temperature value of the London 
penetration depth. Recent measurement of Basov et al. 
§ showed for pure YBCO that A a (0) ~ 1600A and 
A b (0) ~ IOOOA. In the pure limit, the square of the 
inverse penetration depth is given by 
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and is associated with the superfluid density. Here e is 
electron charge, c the velocity of light, is fJ- compo- 
nent of the Fermi velocity, and / the Fermi-Dirac dis- 
tribution function. The electron energy in the normal 
state is £k and E k is the quasiparticle energy in the su- 
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FIG. 2. Solid line is the ratio of a^/cn as a function of 
cutoff angle <j> c for the chain contribution. Dashed lines cor- 
respond to a° and a\ individually in arbitrary units. 



The zero- 



perconducting state with E k = \fe\ 
temperature value of A M is determined solely by the sec- 
ond term in (|2|) and depends only on normal-state pa- 
rameters. By considering an incomplete Fermi-surface 
integration for the chain layer, one obtains 



Ag(0) 
A2(0) 



(vpc = (cos 2 4>)c 



(3) 



Thus, similar to Eq. (|I|) for the resistivity ratio, the intro- 
duction of a cutoff angle will naturally lead to the right 
ratio for the observed zero-temperature anisotropy of the 
penetration depths. (One must add in a plane contribu- 
tion when comparison with experiment is made.) 

Superconducting- State Transport. — Using the Kubo 
formalism, the zero-temperature microwave conductivity 
along the /i-direction for the chain layer is given by 
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where c denotes an incomplete Fermi surface average, 
7 = iY,{)(uj = 0) is the impurity-induced scattering rate, 
and Ak = Ak + Si (a; = 0) is the impurity renormal- 
ized gap at zero frequency with E, the self-energies (see 
later). It is known that, for a gap with d x 2 __ y i -symmetry 
and a complete Fermi surface average, there will be no 
re-normalization effect on the gap (Si = 0), but this is no 
longer the case for a chain. After carrying out the energy 
integration, Eq. fl) is reduced to 



<7 M (0) = iV(0)e 2 



7 2 W7 2 
[ 7 2 + A 2 (0)]l 



(5) 



where AT(0) = m/(2-iTh ) is the electronic density of states 
at the Fermi surface and m is the electron mass for the 
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band. For the a-axis, 7 a (</>) = cos0 and for the 6-axis, 
7b (0) = sin0. 

To show the physics in a transparent way, we first 
ignore the gap renormalization, A(0) — > A(</>) = 
A cos(2</>), and manipulate Eq. (0) to obtain cr a ,b(0) = 
(ne 2 /7rmAo)/ ,b(^c, 7), where n = kp/2ir is the layer 
carrier density and 
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Here, 7 = 7/A0 is the normalized scattering rate. The 
effect of impurities on the renormalization of the gap and 
on the effective scattering rate will be discussed later. In 
Fig. |^, we plot I a ,b for two different choices of 7, namely 
7 = 0.1 and 0.01. For small 7, the function I a is in- 
distinguishable from It, at all values of 4> c . Only when 
7 is large, does Ib starts to deviate from I a and, due 
to the fact that the chain Fermi surface mainly develops 
in the (/) = ±n/2 region and has no contribution from 
the <j) = 0,7r region, <Jb(0) is slightly larger than a a (0). 
More importantly, / Q ,b almost vanish as long as 4> c > 4>q 
(nodal angle) in the small 7 case. (In our present case 
for simplicity, we have 4>q — 7r/4.) 
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FIG. 3. Values of I a and h defined in Eq. (g) in terms of the 
cutoff angle <f) c . The angle (j>o corresponds to the node. Top 
frame is for 7/A0 = 0.1 and bottom frame is for 7/A0 = 0.01. 

The small saturated low-temperature microwave con- 
ductivity in a d a; 2_ i( 2-wave superconductor is recovered 
from Fig. ^ in the case of <p c — ► 0. As is clearly seen 
in Fig. ||, both I a and If, approach one and the con- 
ductivity is saturated to the value a a {0) = o"b(0) 



Ug = ne /(7rmAo), independent of 7, 



of the im- 



purity concentration. It is worth emphasizing that both 
a a (0) and <7b(0) are close to the saturated value <j° s as 
long as 4> c < <f>Q. The latter case means that the gap 
nodes cross pieces of the Fermi surface. In the oppo- 
site case 4> c > $0, the residual conductivity cr a (0) and 
<7b(0) (or I a and lb) are both strongly suppressed. In 



this case, no gap nodes cross the one-dimensional chain 
Fermi surface. This is an effect due solely to the geom- 
etry of the chain Fermi surface and arises even though 
the gap used in our model calculation has pure d x 2_ y 2 
symmetry. Moreover, when <j) c > 4>q and a minimum 
gap A m i n = Arj| cos(20 c )| is seen on the chain Fermi sur- 
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face, one can show that I a < 7r(cos 2 0)c7 2 / (7 2 
and I b < 7r(sin 2 0) c 7 2 /(7 2 + A min ) 3 / 2 . Here A min = 
A m i n /Ao. This means that both I a and Ib are negligibly 
small when 7 <§; A m ; n . 

Thus, for the case of YBCO under consideration here, 
there is a possibility that the gap nodes do not cross 
the entire Fermi surface for the highly one-dimensional 
chain band (4> c > 4>o) and, as a consequence, no satura- 
tion for o~ a ,b(0) occurs in the chain layer. In addition, for 
the same reason that the residual conductivity is strongly 
suppressed in chain layer (for both a and 6-axes), there is 
no significant anisotropy in the residual microwave con- 
ductivity. 

In order to understand what the temperature regime 
is, in which the universal value is valid for YBCO, one 
recalls that the low-temperature microwave conductiv- 
ity, cr(T) ~ cr°(l + T 2 /7 2 ), as given by Hirschfeld et 
al. pi for a g?^^ 2 -wave superconductor with a com- 
plete circular Fermi-surface. Therefore, for the Cu02 
planes only, one needs to have T 7 to achieve uni- 
versality. For the chain layer in which the gap nodes do 
not cross the chain-band Fermi surface, a minimum gap 
A m j n is seen and a result similar to the s-wave case holds 
for the low temperature residual conductivity, namely 
ct(T) ~ a n (Amin / T)e~ A ">j"/ T is predicted (see, for exam- 
ple, Ref. |15|). The chain residual conductivity is thus 
strongly suppressed until the temperature T > A m i n . 
The above picture corresponds exactly to the strong sup- 
pression of tr a (0) and Ob(0) for chain layer discussed be- 
fore. Hence, the chain layer will have little effect on the 
temperature limit where universality holds if A m ; n asso- 
ciated with the chain is larger than 7 associated with the 
plane. 

Now we discuss the effect of impurity scattering on the 
zero- frequency effective scattering rate [7 = iSo(w = 0)] 
and gap renormalization [<5Ak = Ak — Ak = Si (a; = 0)]. 
In the usual t-matrix approximation, the self-energies 
£q(w) = riiGo/ (c 2 - Gq + G\) and Ei(w) = -niG x /{<?- 



Gq + G\) with rii the impurity concentration, c 2 = 1/v'f 
(vi is the impurity potential) , and the integrated Green's 
function Gi(u) = J^ k Gj(k, w). We note in general that 
the impurity renormalization effect should also be in- 
cluded in the band structure <5ek = S3 (in a general sense, 
4> c will also be renormalized in our model). However, this 
is safely neglected. 

In the unitary resonant scattering limit (c = 0), we 
obtain the coupled equations (<5 Ak = 6) 
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FIG. 4. The impurity-induced effective scattering rate 7 
and gap renormalization 8 scaled to Ao (gap maximum in the 
absence of impurity). Solid line is for m/AN(0)Ao = 0.1 and 
dashed line is for n.j/4iV(0)Ao = 0.01 with m the impurity 
concentration. 

where 



\ + (cos 20 + 8) 2 1 

with 8 = 6/A0. In Fig. we show our self-consistent 
solutions for 7 and 8 as a function of 4> c given by Eqs. (|^) 
and (|§|) . We have used two different impurity concentra- 
tions, namely rij/4iV(0)Ao = 0.1 and 0.01. Two features 
are noted in Fig. |[ (i) 7 is strongly suppressed in the 
large 4> c case and (ii) the gap renormalization 8 is al- 
ways negative. The strong suppression of 7 in the larger 
(f> c case is a natural consequence of the evolution of a d- 
wave superconductor towards an s-wave superconductor 
(i.e., a minimum gap A m j n develops in the chain). On the 
other hand, the result that the gap renormalization 8 is 
negative is because the highly ID Fermi surface samples 
mainly the negative portion of the e?-wave gap. 

Analytical results are available in three different 
regimes. In the (f> c — > limit, one recovers the results for 
the complete Fermi surface case (appropriate for Cu02 
plane) that is 7 - (Ao/r) 1 / 2 [r^ 1 ~ rn/N(0) is the 
normal-state scattering rate] and 8 = (i.e., no gap 
renormalization) . When <f) c is small but finite (in the case 
that Fermi surface still sees nodes but the band struc- 
ture has evolved from the tetragonal to the orthorhom- 
bic limit), we find Iq — —7 In 7 and I\ ~ —l/<f> c . Conse- 
quently, 7 ~ A exp(— tA o 2 ) and 8 ~ — 1/(t</> c ). Thus, 
7 is decaying exponentially as <f> c increases and \8\ is pro- 
portional to l/(f> c . When <p c > </>o such that a minimum 
gap A m j n develops on the Fermi surface, we find 7 is small 
(similar to the s-wave case), while 8 ~ —1/[t(tt/2 — 4> c )\- 



Therefore, \8\ increases as <fi c increases and diverges when 
cj) c approaches tt/2. Of course, the latter case is an un- 
physical regime. The three different regimes for 8 are 
seen most clearly in Fig. [|. 

Since 8 < 0, the gap nodes will be effectively shifted 
downward (or equivalently (f>o < 7r/4). As a result, this 
effect increases the chance that the gap nodes will not 
cross the entire chain Fermi surface. In addition to the 
effect of impurity, the appearance of a negative s-wave 
admixture in the d x 2_ y 2-w&ve gap ]l^ , [l4| can naturally 
distort the gap nodes to be less than 45 degree and thus 
leads to a similar consequence. A finite minimum gap on 
the chain Fermi surface is therefore seen to be a robust 
feature of the one dimensional nature of the bands and 
of a d-w&ve gap. Recently, a detailed survey of the ef- 
fect of Ni impurity (substitution for Cu) on YBCO has 
been carried out by Homes et al. fl6f| using optical mea- 
surement. They concluded that the deposition of Ni is 
mainly situated in the chain layer. Hence the change 
with Ni impurity concentration can be used to test our 
model. 
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